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Abstract. The stability and conservation properties of a recently proposed poly- 
merization model are studied. The achiral (racemic) solution is linearly unstable 
once the relevant control parameter (here the fidelity of the catalyst) exceeds a 
critical value. The growth rate is calculated for different fidelity parameters and 
cross-inhibition rates. A chirality parameter is defined and shown to be conserved 
by the nonlinear terms of the model. Finally, a truncated version of the model is 
used to derive a set of two ordinary differential equations and it is argued that these 
equations are more realistic than those used in earlier models of that form. 
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1. Introduction 

The chirality of molecules in living organisms must have been fixed 
at an early stage in the development of life. All life that we know is 
based on RNA and DNA molecules with dextrarotatory sugars. There 
is growing evidence that the RNA world (Woese, 1967; Crick, 1968; 
Orgel, 1968; see also Wattis &: Coveney 1999) must have been preceded 
by a simpler pre-RNA world made up of achiral constituents (Bada, 
1995, Nelson et al, 2000). An alternative carrier of genetic code are 
peptide nucleic acids or PNA (Nielsen, 1993). These can be rather 
simple and are currently discussed in connection with the idea to build 
artificial life (Rasmussen et al., 2003). Furthermore, although PNA can 
still be chiral (Tedeschi et al., 2002), there are also forms of PNA that 
are achiral (Pooga et al., 2001), suggesting that chirality may have 
developed later when the first RNA molecules formed. 

In current proposals to build artificial life, chirality does not seem to 
be crucial. The PNA molecules is proposed to act primarily as charge 
carrier, i.e. a very primitive functionality compared to the genetic code 
in contemporary cells (Rasmussen et al., 2003). At this stage, homochi- 
rality may have been introduced by chance. This is also supported 
by the fact that chiral polymers of the same chirality tend to have a 
more stable structure (Pogodina et al., 2001) and would therefore be 
genetically preferred. 
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Since the introduction of chiral molecules is assumed to take place 
at a stage when there is already growth and self-replication, it is also 
plausible to assume that the existence of chiral molecules has an auto- 
catalytic effect in producing new chiral molecules of the same chirality 
(Kondepudi et al, 1990). This is the basis of the recently proposed 
polymerization model of Sandars (2003); see also Wattis & Coveney 
(2004). The purpose of the present paper is to reconsider this model 
(or a slightly modified version of it) and to analyze its stability behavior 
and conservation properties. We also discuss and illustrate some of the 
salient features of the model in more detail. The model is then compared 
with earlier models of homochirality where the detailed polymerization 
process is ignored and the dynamics of single variables representing left 
and right handed polymers are modeled instead (Frank, 1953; Kon- 
depudi and Nelson, 1984; Goldanskii and Kuzmin 1989; Avetisov and 
Goldanskii, 1993; Saito and Hyuga, 2004). 

In order to appreciate the nature of the many terms in the model of 
Sandars we begin by discussing first the basic principle of the model in 
connection with homochiral polymer growth and then turn to the full 
set of reactions that are included in the model. 



2. Homochiral polymer growth 

In this section we discuss the growth of polymers by adding monomers 
of the same chirality, i.e. we ignore reactions with monomers of the 
opposite chirality. This is conceptually the simplest case, but its equilib- 
rium solution also corresponds to a solution of the full system discussed 
below. We write down the equations for left-handed polymers, but the 
same applies also to right-handed polymers. 

A left-handed polymer of length n is assumed to react with a left- 
handed monomer via the reaction 

Ln + Li^Ln+l (n>l). (1) 

The reaction rate is ks, but since L„ can bind to Li on either side, 
the total reaction rate is 2ks and proportional to the product of the 
concentrations of the two constituents. We denote the concentration 
of Ln chains by so in a volume V the number of L„ chains is 

N.,i = [Ln]V. For n > 3 the number of possible pairs of Ln-i and 
Li is Nn X Ni. A special situation arises for n = 2, because then Li 
is interacting with another Li, and the number of possible pairs is 
only ^Ni{Ni — 1) ~ ^N^. [This problem is familiar from the physics 
of nuclear reactions; see, e.g., Kippenhahn and Weigert (1990).] We 
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therefore introduce the factor an^"^^ defined by 

t III (^) 

(Later we shall use this factor also with a = instead of 1/2.) The 
corresponding contribution to the evolution of the concentration of L„ 
is therefore 

^[^ = ... + 2fc54^/^)[L„]M, (3) 

where the dots denote the presence of other terms that will be discussed 
later. 

Obviously, the concentrations of L„ and Li have to decrease at the 
same rate by the same amount, so 

^ = ...-2ksa^y^\L^][L,l (4) 

^ = ...-2fc5a(V2)[L„_,][Li]. (5) 

In the following we regard n as a general index with 2 < n < N, 
the evolution of [L„] is governed by the difference of two terms (gain 
from Ln-i chains and loss in favor of producing Ln+i chains). The 
production of each L„ contributes to a loss of Li monomers, so the 
right hand side of Eq. (5) becomes a sum over all n. The full set of 
equations is then 



d[Li 
dt 



N-l 

Ql-\l[Li], where Xl = 2ks J^i^n], (6) 



(7) 



^ = 2ks[L,](aW'^[L^.,]-[L^ 

where Ql denotes the production of new Li monomers (see below). A 
corresponding set of equations applies also to right-handed polymers, 
i.e. Ri and Note that Eqs. (6) and (7) obey the conservation law 

dEL 



dt 

where 



QL-2ks[Li][LN], (8) 



AT 

El = J2 ^[Ln] (9) 

n=l 

is the total number of left-handed building blocks. This number reaches 
an equilibrium if the supply of new left-handed monomers, Ql, bal- 
ances the loss associated with reactions involving the longest polymers 
possible for a given value of A'". 
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Equation (7) shows that in the steady state we have = ^[Li] 
for all n > 2. Using Eq. (6), we find Al = ksN[Li], and therefore 

2[Ln] = [Li] = \jQL/ksN (steady state) (10) 

is a possible equilibrium solution. 

New left and right handed monomers are assumed to be continuously 
reproduced from an achiral (racemic) substrate. The rates of regener- 
ation, Ql and Qr, depend on the concentration of the substrate, [S], 
and in some fashion on the relative concentrations of right and left 
handed polymers. So, in general, we write 

Ql = kc[S]{^il + /)Cl + 1(1 - f)CR + Col}, (11) 

Qr = kc[S]{^{l + f)CR + \{l - f)CL + Coi?}, (12) 

where Cl and Cr are some measures of the catalytic effect of the 
already existing right and left handed polymers, and the terms Cql 
and Cor allow for the possibility of non-catalytic production of left 
and right handed monomers - possibly at different rates. (Unless noted 
otherwise, we keep Col = Cqr = 0.) 

The concentration of the substrate is assumed to be maintained by 
a source Q, so we have 

'^ = Q-{Ql + Qr), (13) 

where Ql + Qr = kc[S]{CL + CR + CQL + CoR)- see Eqs. (11) and (12). 
In general, wc expect Cl and Cr to be some function of L„ and 
respectively. Sandars (2003) assumed Cl = [Ln\ and Cr = [-Rjv], 
i.e. the catalytic effect depends on the concentrations of the longest 
possible chains of left and right handed polymers. This assumption 
imposes a dependence on the cutoff value A^, a dependence that should 
preferably be avoided in numerical or other technical considerations. 
The model should for example be stable and consistent in the limit 
when N is infinite. Another option would be to assume Cl = [Lm] 
and Cr = [Rm]-, where M < is a fixed value that is independent 
of the maximum chain length. Both alternatives have the disadvantage 
that [Li] and [i?i] can never grow unless [Lm] or [Rm] are initially also 
finite. While it is plausible that long chains carry more catalytic weight 
than shorter ones, the dependence of the results on the particular choice 
of M seems artificial. (The allowance of finite values of Col and Cqr 
would remove this problem, although in practice both of these values 
should still be quite small.) 
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On may expect that the catalytic properties of the existing left and 
right handed polymers depend on the length of the polymer. The exact 
functional expression for this dependence is not known. It is therefore 
important that a model that explains homochirality is not sensitive to 
the details of the catalytic properties and hence the functional form 
of and Cr. It turns out that the qualitative behavior of the model 
of Sandars is indeed robust in this respect, e.g. a pitchfork bifurcation 
exists in both Sandars' original and in our model. To avoid artificial de- 
pendence on the maximal chain length we chose to let the catalytic 
functions have the following form 

Cl = El, Cr = Er, (14) 

where E^ is given by Eq. (9), and Er is defined analogously. This 
is similar to the choice of Wattis Sz Coveney (2004) who assumed, 
independently of us, Cl = El — [Li] and Cr = Er — [Ri]- 

We now comment on another aspect of the model of Sandars. He 
assumed that in the evolution of [Ln] the loss is not 2ks[Li][LN], as it 
would be if Eq. (7) were applied ton = N, but he introduced an explicit 
linear damping term instead. This implies that the model behaves dis- 
continuously at the end of the chain. We feel that an "extrapolating" 
(continuous) behavior is more reasonable, so we choose to apply Eq. (7) 
also at n = N. 

It is interesting to note that in the continuous limit, Eq. (7) becomes 

which describes waves traveling toward larger n. This is shown in 
Figure 1, where we have perturbed the equilibrium solution (10) by a 
gaussian and have solved Eqs. (6)-(7) numerically. The wave is damped 
and has a speed that is proportional to (ksQ)^/'^, because for the 
steady state background solution [Li] (ks/Q)^^"^ ■ Note that the 
extrapolating boundary condition at n = allows the wave to escape 
freely. 

In this paper we do not adopt the nondimensionalization of Sandars. 
Instead, we note that there are only two physical dimensions in this 
problem: time and volume. Characteristic quantities with the dimen- 
sions of time and volume are (ksQ)-^/^ and {ks/Qy/^, respectively. We 
therefore present all results by explicitly quoting these dimensions. In 
practice this means that from now on we use Q = ks = 1 ss numerical 
values, but we keep the symbols in the equations for clarity. Throughout 
this paper we also assume kc/kg = 1; calculations with different values 
do not seem to affect our results in any important way. 
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Figure 1. Wave- like propagation of a finite amplitude perturbation. The initial pro- 
file is a gaussian. Note the undisturbed propagation of the wave out of the chain at 
n = N. The time difference between the different curves is 20/ [ksQy^^. We have 
shown the first and last times as solid and dashed lines, and all other times as dotted 
lines. The parameters are N = 50 and kc/ks = 1- 

The fact that the equihbrium solution is constant for all n > 2 
implies that this value will decrease for longer choices of N. In that 
sense the solution is never converged. This situation changes when 
we allow the ends of the left-handed polymers to be spoiled by right- 
handed monomers, as done by Sandars (2003). This will be discussed 
in the next section. 



Already 20 years ago, Joyce et al., (1984) showed in an important 
paper describing experiments with template-directed polymerization 
that, once a monomer of the opposite chirality is bound to one end 
of the chain, the polymerization terminates on that end of the chain. 
Sandars (2003) incorporated this effect in his model and showed that 
this can lead to a bifurcation into two possible solutions of opposite 
chirality and hence to homochirality. 



3. Enantiomeric cross-inhibition 
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The full set of reactions included in his model is (for n > 2) 

(16) 
(17) 
(18) 
(19) 

and for all four equations wc have the complementary reactions ob- 
tained by exchanging L R. Following Sandars (2003), we have 
introduced the new parameter kj, which quantifies that rate of enan- 
tiomeric cross-inhibition. The special case A;/ = corresponds to the 
case discussed in the previous section. 

The most important effect of enantiomeric cross-inhibition is that 
a certain fraction of chains becomes spoiled by producing LnRi and 
RnLi polymers. Equation (7) and its complementary equation for right- 
handed polymers suffer therefore a loss proportional to 2ki, so we have 
instead 

^ = 2ks[Li](a^y^)[Ln-i] - [Ln]) - 2ki[Ln][Ril (20) 

^ = 2ks[Ri]{a^^'^\Rn-i] - [Rn]) - 2ki[Rn][Li]. (21) 

These equations allow us to see what happens in the racemic case with 
[R'n\ = [L'n\- In a steady state we have (for n > 2) 

[Ln] = \a~^"-~^\Li] (racemic solution), (22) 

where we have defined a = 1 + fc/Z/cg. In particular, if kj = ks, then 
[Ln] = 2~"'[Li], i.e. [Ln] drops by a factor of 2 from one n to the next, 
except for n = 1 to 2, where it drops by a factor of 4. We should note, 
however, that this solution can be unstable (see Section 4). 

So far, we have not yet considered the evolution equations for the 
concentrations of the mixed terms, LnRi and RnLi. Following Sandars 
(2003), we abbreviate the corresponding concentrations by [Ln-R] and 
[RnL], respectively, i.e. without the subscript 1 on the terminating 
end of the chain. The effect of generating these terms was already 
manifested in Eqs. (20) and (21) through the appearance of the last 
term proportional to 2ki. Nevertheless, we do need to solve for 
and [RnL] explicitly, because the reactions (18) and (19) consume Li 
and Ri monomers, respectively. The evolution equations for [Li] and 
[Ri] are therefore given by 

^ = Ql-Xl[Li], ^ = Q«-A«[i?i], (23) 
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where 

N-l N N-l N 

Al = 2ks Yl [Ln] + 2ki [Rn] + kg ^ [LnR] + kiYl [RnL] , (24) 

n=l n=l n=2 n=2 

N-l N N-l N 

\r = 2ks Y i^n] + 2ki Y [Ln] + % E i^nL] + fc/ ^ [L„i?], (25) 

n=l n=l n=2 n=2 

are the decay rates that quantify the losses associated with the reactions 
(16)-(19), respectively. 

In Eqs. (24) and (25) the concentrations [L^R] and [RnL] enter, so 
we have to solve their corresponding evolution equations (for n > 2) 

= ks[L,] (aW [L^_,R] - [LnR]) +ki[Ri] (2[L„] - [L„i?]) , (26) 

= ks[Ri]{al^^[Rn-iL]-[R^L])+kj[L,]{2[Rn]-[RnL]), (27) 

where the an^ factor turns off the first term for n = 2; see Eq. (2). In 
Eqs. (26) and (27) the first two terms proportional to A;^ correspond 
to the homochiral growth on the unspoiled end, i.e. reaction (18). The 
third term comes from reaction (17) and the fourth term comes from 
reaction (19) and enters here and also in Eqs. (24) and (25) as a loss 
term. For completeness, we note that the corresponding gain enters in 
the evolution equations 

^1^^ = ki[R,][L^R], &^ = ki[L,][RnL], (28) 

which arc not explicitly required for constructing a solution, because 
these polymers no longer react with the monomers. Nevertheless, solv- 
ing Eq. (28) simultaneously with Eqs. (20) and (21) and Eqs. (23)-(27) 
can be useful for monitoring the evolution of the net chirality; see 
Section 5. 

Note that, in contrast to the equations given by Sandars (2003), 
the truncation levels for the terms [L„], [L„i?], and [RLnR] are here 
the same, i.e. n < N for both terms, whereas in the model of Sandars 
the longest chain has n = N — 1, and the longest i?iL„i?i chain 

has only n = N — 2. The reason we need to keep the same truncation 
levels for all three types of polymers is that we want to ensure that 
the behavior near the end of the chain (n = N) does not deviate from 
the behavior elsewhere (n < A''); see the discussion in Section 2. For 
example, to ensure continuous behavior of [Ln] at n = N we need to 
keep the term — 2/c/[Ljv][-Ri] in Eq. (20). This term, however, is the 
loss resulting from the gain of [LnR] , so we have to keep the evolution 
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Figure 2. [Ln] (left) and [LnR\ (right) of equilibrium solutions for different values 
of /. For / = 1 we have [LnFt\ = 0, which cannot be seen in the logarithmic 
representation. 



equation for this term as well. Furthermore, the evolution equation for 
this term involves, in turn, the term —ki[Ri][Li^R], which is the loss 
corresponding to the gain of [RLj\jR\. If one regards the truncation 
level N as an unrealistic feature of the model, as we do, then all three 
polymer types should be truncated at the same level. 

In Figure 2 we show [L„] and [i^^i?] for a number of equilibrium 
solutions for different values of / and kj/ks = 1- The corresponding 
values of [Rn] and [RnL] are small and not shown, except when / = in 
which case the solution is fully racemic with = [L„] and [RnL] = 
[LnR] and is simply 

[LnR] = (n - 1)[L„] = {n- 1)2"" [Li]. (29) 

For / = 1 the solution is given by Eq. (10). 

For ki/ks = 0.1 the results are similar to those for kj/ks = 1 
provided / > 0.8. For / < 0.7, however, the solution is fully racemic 
and therefore the curves are independent of /. This racemic solution is 
similar to the case kj/ks = I and / = 0.8 that is shown in Figure 2. 



4. Stability of the racemic equilibrium 

A realistic model of homochirality must also have an achiral (racemic) 
equilibrium solution. It is generally anticipated that this racemic solu- 
tion can be destabilized in the presence of catalytic reactions (Frank, 
1953; Avetisov and Goldanskii, 1993). If the probabilities of adding left 
and right handed monomers to a homochiral polymer are equal, i.e. 
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Figure 3. Growth rate (left) and bifurcation diagram showing a classical pitchfork 
bifurcation (right) as a function of fidelity for kj/ks = 1, and N = 50. The dashed 
line indicates an unstable solution. 



if kj = ks, the racemic solution given by Eq. (29) is also a possible 
solution for other values of the fidelity than / = 0, but it may of course 
be unstable. 

We have carried out a numerical stability analysis by adding a small 
(10~^) relative perturbation to the value of [Li] of the racemic solution. 
It turns out that for certain values of the fidelity / the departure of 
[Li] from the racemic equilibrium solution, (^[Li], growth exponentially 
in time like e^*. In Figs 3 and 4 we plot A obtained from the slope of 
the graph of ln(5[-Li](t) during the exponential growth phase, i.e. before 
a new nonlinear equilibrium is attained. In Figs 3 and 4 we also plot 
the corresponding chiral polarization parameter, ij, as a function of /. 
Here we have chosen to define rj as 

r] = iER- El) / (Er + El) . (30) 

It turns out that for kj/ks = 1 the racemic solution is unstable when 
/ > 0.39, and for kj/ks = 0.1 it is unstable when / > 0.735. The 
transition from an achiral to a chiral solution is a typical example of 
a pitchfork bifurcation; see Figs 3 and 4. This result is in qualitative 
agreement with Sandars (2003) who found that for kj/ks = 1 the 
critical value of / is around 0.21. The differences in the numerical values 
are explained by differences in the model (e.g., the coupling to the 
substrate and the length of the maximum polymer length). 

The growth rate of the instability is important for determining the 
time it takes for an almost racemic solution to become homochiral (or 
at least non-raccmic for / 7^ 1). When kj/ks = 1, the growth rate A is 
around 0.5, but it becomes significantly smaller when the value of kj 
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Figure 4- Same as Figure 3, but for ki/ks = 0.1. 

is reduced. This shows expHcitly that homochirahty emerges as being 
due to enantiomeric cross-inhibition. 



5. Conservation of chirality 

For homochiral growth the relevant conservation law is given by Eq. (8) 
for El, and similarly for Efi. In general, however, because of the in- 
teraction with left and right handed monomers, there are no longer 
separate conservation laws for El and Instead, the complete set of 
equations, Eqs. (20) and (21) together with Eqs. (23)-(25), satisfies 



—AE = AQ- AA, 
dt 



(31) 



where AQ = Qr — Ql and AA = Ar — Al are the net input and output 
rates of chirality, respectively, and AE = Er — El is the total chirality, 
where 



N 



N 



N 



Er = J2 n[Rn] + - ^)[RnL] + - '^)[LRnL], (32) 

2)[RLnR], (33) 



n=l 
N 



n=2 
N 



n=3 
N 



n=l n=2 n=3 

denote the total numbers of right and left handed building blocks 
(or enantiomers) , where opposite enantiomers are counted such that 
they annihilate enantiomers of the opposite chirality. The loss terms 
resulting from the finite truncation level, N, arc denoted by 



Ar = 2ksN[Ri][RN] + ksiN- 1)[Ri][RnL], 



(34) 
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Al = 2ksN[Li][LN] + ks{N - 1)[Li][LnR]. (35) 

In order to evaluate the quantities Er and El we have to integrate 
the evolution equations (28) for the production of terminally spoiled 
polymers - even though they undergo no further evolution. In a sense 
the integration of the terminally spoiled polymers acts only as counters 
that keeps track of the number of polymers that are lost during the 
polymerization process. 

The expressions for Ef{ and El involve sums over [LRnL] and [RLnR], 
but since these quantities do not occur on the right hand sides of the 
governing evolution equations, their values are not constrained by the 
dynamics and depend on the initial conditions and continue to evolve 
in time even though the system may have reached an equilibrium. The 
so defined net chirality can therefore not be used to characterize a 
particular solution, and we have to restrict ourselves either to Er and 
El, or to Er and El, which are defined by taking only the first two 
sums in Eqs. (36) and (37), i.e. 

N N 

Er=Y1 ^i^n] + - ^)[RnL], (36) 

n=l n=2 

TV TV 

El = Y1 ^i^n] + E(" - l)[^ni?]- (37) 
n=l n=2 

In Table I we list the resulting values of r/, defined in Eq. (30), an 

analogously defined fi = {Er - El)/{Er + El), AE = Er - El, 
and AE = Er — El- We also give the mean polymer lengths, Nr = 
J2n[Rn]/J2[Rn] and Nl = J2n[Ln]/ J2[Ln], of right and left handed 
polymers. 



6. Comparison with other models 

The polymerization model of Sandars (2003) is significantly different 
from all the previously proposed models of homochirality that ignore 
the detailed polymerization process by only describing some scalar 
quantities, say x and y, that are representative of the number of left and 
right handed polymers. In the papers by Saito and Hyuga (2004) it was 
shown that neither linear nor nonlinear autocatalytic behavior suffice to 
produce homochirality, and that a backreaction term is needed. Their 
model equations are 

X = x\\ - r) - ex ^ 

y = y\l-r)-ey ^ ' 
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Table I. Numerical results for rj, tj, AE, and AE for different values of /. The ± 

indicates that these values can have either sign. The last column gives the typical 
value of TV necessary for obtaining a converged result representing the limit N — > oo. 
For / = 1 the results for AE {— AE) and Nr do not converge and we give the 
analytic expression for arbitrary A'^ instead. 
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where r = x + y and e is the feedback parameter. For e = there 
is a continuous range of solutions along the line r = x -\- y = I, i.e. 
homochirality docs not emerge unless the initial condition is already 
homochiral. For finite (but small) values of e there are two nontrivial 
stable fixed points. (The trivial solution, x = y = 0, is always a stable 
fixed point in this model.) 

The model of Saito and Hyuga (2004), hereafter the SH model, does 
capture the expected behavior, but it remains unsatisfactory in that its 
functional form has been introduced ad hoc. It is therefore desirable to 
derive simple model equations based on the polymerization equations of 
Sandars (2003). It turns out that, without changing the basic properties 
of the model, a minimal version is still meaningful for N = 2, and that 
the equations for the semi-spoiled polymers, [L2i?] and [R2L], can be 
ignored (as already done by Sandars). Thus, we only solve Eqs. (20) 
and (21) together with Eqs. (23)~(25). Following Sandars (2003), we 
also assume that Cl = [L2] and Cr = [R2] (instead of Cl = El and 
Cr = Er, which would yield more complicated expressions). A further 
simplification can be made by regarding [L2] as a rapidly adjusting 
variable that is enslaved to [Li] (and similarly for [i?2])- This technique 
is also known as the adiabatic elimination of rapidly adjusting variables 
(e.g., Haken, 1983). Equation (20) becomes 

= ks[Lif - 2[L2] {ks[Li] + ki[Ri]) , (39) 
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/=0.95 /=0.501 




Figure 5. Phase diagram showing the trajectories of solutions of Eq. (42) for two 
different values of /. The starting points of each trajectory are marked by small dots 
and stable fixed points axe marked by big dots. 



which is solved for [L2] (and similarly for [R2]), which in turn couples 
back to the equations for [Li] and [Ri] via Qi and Qji. Finally, we 
also treat the substrate [S] as a rapidly adjusting variable, i.e. we have 
kc[S] = Q/{[L2] + [R2\)- We emphasize that the adiabatic elimination 
does not affect the accuracy of steady solutions. It is convenient to 
introduce new dimensionless variables, 

x=[Ri]{2ks/Q)'/'', y=[Lr]{2kslQf'\ r = t{Qks/2)^/'' . (40) 

In order to compare first with the SH model we restrict ourselves to the 
special case kj/ks = f = 1, which leads to the revised model equations 

X = x^r^ - rx, 
y = y /r^ -ry, 

where dots denote derivatives with respect to r and r = x + y and 
= + y'^ has been introduced for brevity. Equations (41) resemble 
the equations of the SH model in that both have a quadratic term 
proportional to x^ (or y^), which is quenched either by a 1 — r factor 
(in the SH model) or by a 1/f^ factor in our model. Furthermore, both 
models have a backreaction term proportional to —x (or —y), but the 
coefficient in front of this term (e in the SH model) is not constant but 
equal to r. 

In the general case with A;//A;5 7^ 1, / 7^ 1, as well as finite values of 
Cqx = CoR{2ks/Qy^'^, and Coy = CoL{2ks/Q)^^'^, the equations read 

X = {px^ + qy^ + Cox)/r^ - r^x, 
y = ij>f + qx^ + C'o2/)/r^ - Tyy, 
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Figure 6. Imperfect bifurcation obtained by solving Eq. (42) for Cqx 
Coy = using the Newton-Raphson method. 



0.001 and 



where we have introduced the abbreviations rx = x + ykj/ks, 



P = SP' 



+ y + Cqx + C'oy, 



y + xki/ks, = xy2rx, r = y 
p = (1 + /)/2, and g = (1 - /)/2. 

In Figure 5 we show trajectories of solutions of Eq. (42) for two dif- 
ferent values of / in an (x, y) phase diagram. Note that all equilibrium 
solutions lie on the line r = 1. This property allows us to calculate 
equilibrium solutions for general values of /. Inserting y = 1 — x yields 
a cubic equation of which one solution is always x = 1/2. This reduces 
the problem to a quadratic equation with the solution 



1 



1 



for / > 1/2, 
otherwise. 



(43) 



Linearizing the equations around the racemic solution, x = y = 1/2, 
yields the growth rate 

A = 2/ - 1. (44) 

In agreement with our numerical results for large values of N, this equa- 
tion gives a linear dependence of the growth rate on the fidelity. This 
result also shows that for /< 1/2 perturbations decay exponentially. 

In the presence of a biased, non-catalytic generation of monomers 
(finite Cqx or Coy with Cqx 7^ Coy) there is no longer a perfectly racemic 
equilibrium solution. The sign of rj for the solution for / = depends 
on the sign of Cqx — C'oj/- Along this solution branch rj goes further away 
from zero in a continuous fashion until / = 1. At some value of / a 
pair of new solutions emerges, one is stable and the other one unstable. 
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but both have the opposite sign of r/; see Figure 6. Among these new 
branches, the stable one can only be reached via a finite amplitude 
perturbation. This behavior is called an imperfect bifurcation and has 
long been anticipated in this context (Kondepudi and Nelson, 1983; 
Kondepudi et al, 1986; Goldanskii and Kuzmin, 1989). 

The steady solutions shown in Figure 6 have been obtained by solv- 
ing Eq. (42) using the Newton-Raphson method. This method allows us 
to find both stable and unstable solutions. Near the bifurcation point 
the diagram is extremely sensitive to the addition of a bias parameter. 
It is remarkable that for a value as small as Cqx = 10~^ a relatively 
large gap has been produced in the bifurcation diagram. 

Finite values of Cqx and Coy could result from physical influences, 
for example polarized synchrotron radiation from neutron stars (hut see 
Bonner, 1999), UV radiation in star-forming regions (Bailey, 2001), or 
the parity violation of the electroweak force (e.g., Hegstrom, 1984). In 
all these cases the expected effect is however very small (Bada, 1995). 
We emphasize, however, that the main reason for homochirality is the 
instability of the racemic (or nearly racemic) solution, which is hardly 
modified by a finiteness of Cqx or Coy. 

7. Conclusions 

The origin of homochirality has long been thought to be the result 
of a bifurcation process that can vastly amplify a very small random 
enantiomeric excess which can then prevail forever. Generic model 
equations reproducing the expected bifurcation behavior have so far 
mostly been proposed on an ad hoc basis. It was therefore difficult 
to establish a connection between model and reality. According to the 
work of Saito and Hyuga (2004) one expects two effects to be important: 
nonlinearity and backreaction. However, the functional form of these 
terms remained open. Furthermore, the meaning of non-perfect cat- 
alytic fidelity and enantiomeric cross-inhibition within the framework 
of the model were not clear. In the present paper we have established 
a direct connection between the more detailed polymerization model 
of Sandars (2003) and the simpler model equation approach with only 
two ordinary differential equations. In particular, the present work has 
confirmed that the relevant nonlinearity is indeed quadratic (as in the 
SH model), but it is not quenched like 1 — r, but rather like 1/f^, where 
r and r are measures of the total concentrations of monomers (both 
right and left handed). Furthermore, the feedback coefficient is not a 
small constant, as in the SH model, but it is itself proportional to r. 
More importantly, imperfect fidelity and enantiomeric cross-inhibition, 
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as well as the effects of a weakly biased non-catalytic production of 
new monomers, have a quantitative meaning within the framework of 
the reduced model. 

For a more quantitative comparisons of the polymerization process 
with experiments the full set of equations of Sandars (with the revi- 
sions discussed above) is to be preferred. A number of features that 
can only be captured by the full model. An example is the wave-like 
propagation in the distribution of homochiral polymers. An experi- 
mental confirmation would help to quantify the growth coefficient ks 
characterizing the probability that a polymer grows by a monomer 
of the same chirality. On the other hand, the growth coefficient for 
enantiomeric cross-inhibition, kj, determines primarily the minimum 
fidelity parameter, /, above which bifurcation and hence homochiral 
growth is at all possible. It is indeed quite remarkable, that the main 
reason homochiral growth occurs is that binding with a wrong enan- 
tiomer spoils further polymerization on the corresponding end of the 
chain. This leads to competition which is always a key feature of natural 
selection processes such as these. 

Homochirality in living organisms is a singular phenomenon. Non- 
living chemical systems do in general not have a preferred chirality. In 
the models presented in this paper this is reflected in Figs 3 and 6. 
The region of the phase diagram displaying homochirality is character- 
ized by high fidelity, i.e. high auto-catalytic accuracy. The fidelity is 
expected to be significantly higher in living systems. When an organ- 
ism dies the auto-catalytic polymerization stops and as a consequence 
the fidelity is sharply decreased. The characteristic behavior of the 
polymerization changes from the chiral to the raccmic region of the 
phase diagram. The relaxation of the system from the homochiral to 
the racemic state is often very slow. It was in fact suggested by Hare 
and Mitterer (1967) and later by Bada et al, (1970) that racem- 
ization of amino acids in fossil material could be used as a dating 
method. Unfortunately it has turned out that the rate of racemization 
is strongly temperature dependent, which tends to make this dating 
method unreliable. 
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